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Abstract 

The shape and tails of partial distribution functions (PDF) for 
a financial signal, i.e. the S&PSOO and the turbulent nature of the 
markets are linked through a model encompassing Tsallis nonexten- 
sive statistics and leading to evolution equations of the Langevin and 
Fokker-Planck type. A model originally proposed to describe the in- 
termittent behavior of turbulent flows describes the behavior of nor- 
malized log-returns for such a financial market index, for small and 
large time windows, both for small and large log-returns. These turbu- 
lent market volatility (of normalized log-returns) distributions can be 
sufficiently well fitted with a x^-distribution. The transition between 
the small time scale model of nonextensive, intermittent process and 
the large scale Gaussian extensive homogeneous fluctuation picture is 
found to be at ca. a 200 day time lag. The intermittency exponent (k) 
in the framework of the Kolmogorov log-normal model is found to be 
related to the scaling exponent of the PDF moments, -thereby giving 
weight to the model. The large value of k points to a large number 
of cascades in the turbulent process. The first Kramers- Moyal coeffi- 
cient in the Fokker-Planck equation is almost equal to zero, indicating 



"no restoring force". A comparison is made between normalized log- 
returns and mere price increments. 
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1 Introduction 



The time lag dependent price increments, returns, log-returns, normalized 
log-returns of financial market indices, stocks and foreign currency exchange 
markets are known to be non-gaussian distributions and rather exhibit fat- 
tailed power-law distributions [1, 2, 3, 4, 5]. The origin of the so called large 
volatility characterized by such fat tailed distributions is a key question; the 
fat tails in such data are thought to be caused by some "dynamical process" 
through a hierarchical cascade of short and long-range volatility correlations, 
though Gopikrishnan et al. consider that correlations and tails have different 
origins [5]. Destroying all correlations, e.g. by shuffling the order of the 
fluctuations, is known to cause the fat tails almost to vanish. It is still an 
open question whether both the fat-tailed power-law of partial distribution 
functions (PDF) of the various volatilities and their evolution for different 
time delays in financial markets can be described. 

The fat tails indicate an unexpected high probability of large price changes. 
These extreme events are of utmost importance for risk analysis. They are 
considered to be a set of strong bursts in the energy dissipation of so called 
clusters of high price volatility. In so doing the PDF and the fat tail event 
existence are thought to be similar to the notion of intermittency in turbu- 
lent flows [6]. Indeed employing the Fokker-Planck equation approach [7] 
recent studies [8, 9, 10, 11] have shown that the dynamics of a market results 
from a flow of information between long and short time scales. Since the 
distributions of returns obey a Fokker-Planck equation, the time evolution 
of the price increment signal {Ay) measured for a time lag At is governed 
by a Langevin equation [7] 



the drift D^^^ and diffusion D^'^^ coefficients being those of the Fokker-Planck 
equation [8, 9, 10, 11]. It is often assumed that r]{At) is a correlated noise 
with gaussian statistics. Thus such a dynamics may be analogous to the 
dissipation of energy from large to small spatial scales in three dimensional 
turbulence as pointed out already in [6, 12, 13]. 

On the other hand, the non-Gaussian character of the fully developed 
turbulence [14] has been linked to nonextensive statistical physics [15, 16, 
17, 18, 19, 20, 21, 22, 23]. Whence recently there has been a large number 
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L>(^) {Ay (At) , At) + v{At) ^D(^) {Ay (At) , At) , 
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of studies, e.g. [6, 17, 24, 25, 26, 27] of financial markets employing the 
nonextensive statistics including those involving fully developed turbulence 
approach as in [14]. The nonextensivity, i.e. some anomalous scaling of 
classically extensive properties like the entropy, is linked to a single parameter 
q, e.g. in the Tsallis formulation of nonextensive thermostatistics. 

In this paper on the study of the behavior of a financial index, i.e. the 
S&P500, on large time windows, we apply as in [24, 27] for short time 
windows, a recently suggested model of hydrodynamic turbulence that serves 
as a dynamic foundation for nonextensive statistics [19, 20, 21]. Indeed long 
time lag effects must be also investigated. Furthermore, it is known that 
some distinction must be made concerning the type of financial market signal 
which is examined. We will compare results based on price increment and 
normalized log-return^ time series. 

In Sect. 2, wc describe the distribution of returns for the daily closing 
price signal of the S&P500 index for the time interval between Jan. 01, 1980 
and Dec. 31, 1999, thus a series of =5056 data points. Daily closing 
price values of the S&P500 index for the period of interest were downloaded 
from Yahoo web site ^. We characterize the tail(s) of the distribution for 
various Af s, i.e. from 1 up to 40 days, and will observe the value of the 
PDF tails, for such time lags, outside the best gaussian through the data. In 
Sect. 3, we calculate the power law exponents characterizing the integrated 
distribution of the (normalized log-) returns over different time lags for the 
S&P500 index daily closing price through a detrended fluctuation analysis 
and a power spectral density analysis point of view. Results are compared 
to shuffled data for estimating the value of the error bars. In Sect. 4, Tsallis 
statistical approach is outlined, and distributions of (normalized log-) returns 
for time lags between At = 1 to 40 days are examined. It is found that the 
g- value of the nonextensive entropy converges to a value = 1.22 for At — 
40 days, starting with q — 1.39 for At — 1, values similar to those reported 
for the intraday evolution of other financial indices, e.g. NASDAQ in [28] 
and slightly lower than those for S&P500 minute data [24]. The probability 
density /At(/5) of the volatility (3 in terms of the standard deviation of the 
normalized log returns of the S&P500 for different time lags is found to obey 
the x^-distribution. The intermittency exponent (k) of the Kolmogorov log- 

^ throughout this report the natural logarithm will be always used. 
'^http://finance. yahoo. com 
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normal model is found to be related to the scaling exponent of the PDF 
moments, -thereby giving weight to the model. The large value of k points 
to a large number of cascades in the underlying turbulent process. 

In Sect. 5, the usual Fokker-Planck approach for treating the time- 
dependent probability distribution functions is summarized and coefficients 
governing both the Fokker-Planck equation for the distribution function of 
normalized log returns and the Langevin equation for the time evolution of 
normalized log returns of daily closing price signal of S&P500 are obtained. 
We will notice that there is " no restoring force" . 

Therefore we present for the first time a coherent theory linking the shape 
and tails of partial distribution functions for long and short time lags of the 
daily closing values of a financial signal and connect the often suggested tur- 
bulent nature of the markets to a model encompassing nonextensive statistics 
and evolution equations of the Langevin and Fokker-Planck type. We are 
aware that the number of data points of the time series (A^ =5056) might 
seem quite small with respect to other studies involving millions of data 
points. Some previous work had indicated the possible use of 5000 or so data 
point series in order to obtain scaling arguments and ingredients for models 
[29, 30]. Clearly the relative error bars or confidence interval being roughly 
proportional to A^^^^ have to be taken with caution. Thus the conclusions 
universal value might be debated upon. Nevertheless one positive aspect 
might be that scaling effects are not too sensitive to N. We do warn however 
throughout the report that some caution might also be taken concerning the 
stationarity of the data. These caveat might only be resolved through further 
work. 

2 Distribution of returns 

There are several ways of calculating the returns in a financial market. A 
simple one represents price increment Ay{t, At) or difference between the 
value of the price signal y{t) at time t + At and its value at time t. Log- 
returns, i.e. the logarithm of the price ratio y{t) = ln[y{t + At)/y{t)] 
are also used sometimes. Below we consider the normalized log-returns 
Z{t, At) = {y{t)— < y >At)/c"At, where < y >At denotes the average and ctaj 
the standard deviation of y{t) for a given At. The normalized log-returns 
Z{t, At) depend on the time t and the time lag At. However, in order to sim- 
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plify the notations and whenever possible without leading to confusion and 
misunderstanding we will drop the explicit writing of one or both variables. 
Daily closing price values of the S&P500 index for the period between Jan. 
01, 1980 and Dec. 31, 1999 will serve as a standard financial signal y{t), thus 
N =5056 data points. 

The distribution of the normalized log-returns Z{t, At) of the daily closing 
price signal for S&P500 index for the period between Jan. 01, 1980 and Dec. 
31, 1999, for At — 1 day are plotted in Fig. 1. A fit is first attempted with 
a Gaussian distribution for small values of the increments, i.e. the central 
part of the distribution. This central part distribution is well fitted with such 
a Gaussian type curve within the interval Z G [—3, 3] but departs from the 
Gaussian form outside this interval. The negative and positive tails of the 
distribution outside the Gaussian curve are found both to be equal to -3.1. 
In the case At > 1 day, it is observed that the best Gaussian range is the 
same as for At = 1 day (Fig. 2) but the outside tail values, as estimated for 
the various At's of interest, arc slightly different, and found to decrease with 
At; they are reported in Table 1. These findings (At- independent Gaussian 
range and tail exponent behavior) are at odd with the expectation that the 
PDF tends toward a Gaussian for large At. Some Bayesian-like analysis 
of the PDF's, i.e. allowing for the expected Gaussian width behavior, has 
been done, with the appropriate conclusion. However, the error bar on the 
various widths do not allow for a statistically convincing evidence through 
the comparison of variance classical test. Same for the tail exponents which 
are obtained from a very small number of data points. Therefore it seems 
appropriate to pursue further the PDF analysis through other techniques 
that one allowing to extract a PDF tail from the difference between a raw 
data histogram and a central region Gaussian fit. 

Notice that the Oct 19, 1987 and Oct. 27, 1997 crashes, as studied 
elsewhere [31, 32] are represented by isolated dots at Z — —23 and and 
Z = —7 respectively^. The value of Z = —8.7 represents the aftershock 
crash of the Oct. 26, 1987 [32]. The analysis presented in the present study 
is not designed to capture such extreme events nor their effects. 

^Although the absolute value of the S&P500 drop in price is of order of 60 units in 
both crashes, since the price has increased and due to the nonlinearity of the logarithm the 
value of Z is much smaller at the crash in 1997. 
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3 Time correlations 



There are different estimators for the long and/or short range dependence 
of fluctuations correlations [33]. In all cases it is useful to test the null 
hypothesis debated [4, 5] whether the fat tails are related to or/and caused 
by long-range volatility correlations. Destroying all correlations by shuffling 
the order of the fluctuations, is known to cause the fat tails almost to vanish. 
A Kolmogorov-Smirnov test (not shown) on shuffled data has indicated us 
the statistical validity of the numerical values and the statistically acceptable 
meaning of the displayed error bars. 

Through the (linearly) detrended fluctuation analysis (DFA) method, see 
e.g. [34], we show first that the long range correlations of daily closing price 
signal of S&P500 for the time interval of interest, are brownian-like. The 
method has been used previously to identify whether long range correlations 
exist in non- stationary signals, in many research fields such as e.g. finance 
[29, 30], bioinformatics [35], cardiac dynamics [36] and meteorology [37, 38, 
39] . The DFA concepts are therefore not repeated here. For an extensive hst 
of references see [34]. Briefly, the signal time series y{t) is first integrated, 
to "mimic" a random walk Y{t). The time axis is next divided into k + 1 non- 
overlapping boxes of equal size size n; one looks thereafter for the best (linear) 
trend, z{n), in each box, and calculates the root mean square deviation of 
the (integrated) signal with respect to z{n) in each box. The average of such 
values is taken at flxed box size n in order to obtain 



F{n) 



I 1 (fc+l)n \ 

(- [Y{i)-z{i)f). (2) 



The box size is next varied over the n value. The resulting function is ex- 
pected to behave like ti^+^dfa indicating a scaling law characterized by a 
(Hurst) exponent Hofa- For the (integrated) daily closing price signal of 
the S&P500 index, a scaling exponent 1 + Hdfa = 1-52 ± 0.01 is found (Fig. 

3) in a scaling range extending from about 1 week to about 250 days, i.e. 1 
year. 

Along the same line of thoughts the scaling properties of the normal- 
ized log-returns Z{t,At) = {y{t)— < y >At)/o'At have also been tested for 
different time lag values, i.e. At = 1,5,10,15,20, 25,30,35,40 days (Fig. 

4) . The DFA function, as defined here above, of the integrated normalized 
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log returns for time lag 1 days behaves as white noise and has a Hausdorff 
measure equal to zero (later sec its power spectrum in Fig. 6). However 
non trivial scaling properties occur for the series of normalized log returns 
as soon as At > 1 day. The values of the scaling exponents and the maxi- 
mum box size rix (in days) for which the scaling holds for each DFA-function 
are given in Table 2, while the DFA-functions together with fitting lines are 
plotted in Fig. 4. The values of the Hausdorff measure of the normalized log- 
returns signals varies with At from Hdfa — 0.27 ± 0.04 for At — 5 days to 
Hdfa = 0.43±0.01 for At = 40 days. Recall that Hdfa = 0.5 corresponds to 
Brownian motion. The value of the maximum box size for which the scaling 
holds nx is related to the periodicities of the normalized log-returns signals 
defined by the value of the time lag as Ux ~ 3. 5 At. The data and the power 
law fit of this functional dependence are plotted in the inset of Fig. 4. The 
value of the slope 1.0 is the same as the one found by Hu et. al. [34] when 
studying the effects of sinusoidal trends and noise on the (so called second 
order in [34]) DFA technique. 

The power spectrum of the daily closing price signal of S&P500 S{f) ~ 
f"^ with spectral exponents //i = 2.41 ± 0.06 and 112 — 1-95 ± 0.03 with a 
scale break at 250 days is shown in Fig. 5. The scaling properties of the 
power spectrum of the shuffled daily closing price signal of S&P500, in which 
e.g. the amplitudes are randomly shuffled are shown in the inset of Fig. 5. 
Such a scaling spectral exponent = is the signature of a white noise like 
behavior. Recall that /j, — 2.0 corresponds to Brownian motion. 

We have also checked for scaling behavior and possible periodicities in the 
power spectrum of the time series of the normalized log returns Z{t,At) = 
{y{t)— < y >At)/o'At for different (selected) values of the time lag At = 
1,5,20,35, 40 days (Fig. 6). A white noise like behavior of the power spec- 
trum of such returns always occurs for 1// < At days; e.g., dashed line in 
Fig. 6 for / < 1/128 and At = 40 days. This is in accordance with the 
results of the DFA analysis (Fig. 4 and Table 2). A scaling behavior is found 
at large frequencies / satisfying the relationship /x = 277/)^^ + 1, as indicated 
in Fig. 6, e.g. by the dashed line with slope 11 = 1.86, for / > 1/128 for the 
case At = 40 days. 

Periodicities in the power spectrum of the normalized log return time 
series for At > 1 day were expected to be found since these periods are 
somewhat embedded into the time series by the way they are obtained and 
the Fourier transform technique. It is easily observed that the maxima and 
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the minima of the spectrum correspond to harmonics and subharmonics of 
1/At. 



4 Tsallis statistics 

Based on the scahng properties of multifractals [40], TsaUis [15, 41] proposed 
a generahzed Boltzmann-Gibbs thermo-statistics through the introduction of 
a family of non-extensive entropy functional Sq given by: 

Sg = {l- j P{X, tydX^ , (3) 

with a single parameter q and where k is a normalization constant. The main 
ingredient in Eq.(3) is the time-dependent probability distribution p{x,t) of 
the stochastic variable x. The functional is reduced to the classical extensive 
Boltzmann-Gibbs form in the limit of q — 1. The Tsallis parameter q char- 
acterizes the non-extensivity of the entropy. Subject to certain constraints 
the functional in Eq.(3) seems to yield a probability distribution function of 
the form [6, 15, 19, 24, 27] 



CPo2a{q-l)\x 



2q ~) (q-l) 



2a -{q- 1) 



for the stochastic variable x, where 



r(^)r(^-^) 

in which C is a constant and < a < 1 is the power law exponent of the 
potential U{x) = C|a;p" that provides the "restoring force" F{x) in Beck 
model of turbulence [19, 20, 21, 23]. The latter is described by a Langevin 
equation 

nnr 

- = -^F{x) + R{t) (6) 

where 7 is a parameter and R(t) is a gaussian white noise. A non-zero value 
of 7 corresponds to providing energy to (or draining from) the system by the 
outside [42] The parameter (3o in Eq.(4) and (5) is the mean of the fluctuating 
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volatility f3, i.e. the local standard deviation of over a certain window of 
size m [6]. We will use this model assuming that the normalized log returns 
Z{t,At) represent the stochastic variable x, as in Eq.(6), or Ay in Eq.(l). 
We will search whether Eq.(4) is obeyed for x = Z{t,At), thus studying 

p[x) = p^t{Z) for various time lags At. 

,jTist as in Beck model of turbulence"^ [19, 20, 21] wc assume that the 
volatility P is x^-distributed with degree u (see another formula in [23]): 

where F is the Gamma function, /3o =< f3 > and the number of degrees of 
freedom u can be found from: 

(8) 



< /32 > - < >2- 

The TsaUis parameter q satisfies [19] 

q = l + ^—. 9 
av + 1 

To justify our assumption that the 'local' volatility of the normalized log 
returns Z{t, At) is of the form of x^-distribution, we checked the distribution 
of the normahzed log returns of the daily closing price of S&P500. We have 
calculated the standard deviation of the normalized log returns within various 
non-overlapping windows of size m, ranging from 25 to 1000 days 



m 



\ 



{k+l)m I (fe+l)m 

E ^'w- - E m\ (10) 



i=km+l \ t=km+l 



In doing so we have a various number of Ai non-overlapping windows 
for various time lags At, and have searched for the most efficient size of the 
window in order not to loose data points and therefore information. The 
resulting empirically obtained distributions of the 'local' volatility (Eq.(lO)) 
of normalized log returns for the different different time lags of interest are 
plotted in Fig. 7 for an intermediary case m — 32. The values of the 

^The approach used here was recently suggested to be an appropriate model of hydro- 
dynamic turbulence for financial markets in [27]. 
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degree V of the ^^-distribution are then obtained using Eq. (8). The spread 
[f^mim Pmax] of the local volatihty /? decreases with increasing the time lag as it 
is expected from a x^-distribution function due to the exponential function 
in Eq. (7) for large values of the degree of freedom u. The value of u 
much varies as a function of m and the time lags considered. The fits are 
always excellent. However the (3q and u values are quite dependent on the 
parameters used in the numerical analysis. Based on these results, e.g. Fig. 
7, it can be accepted that the (turbulent market) model /^-distributions can 
be sufficiently well fitted for our purpose with a x^-distribution, thereby 
justifying the initial assumption.^ 

In order to investigate the impact of the a parameter on the tail behavior 
of the Tsallis type distribution function we tested Eq. (4) for fixed q in 
two cases : for a time lag At — 1 day and q — 1.39 (Fig. 8a) and for a 
time lag At = 40 days and q = 1.22 for a = 1,0.9,0.8 (Fig. 8b). Next we 
tested Eq. (4) for fixed a = 1 and varying q : for a time lag At = 1 day 
and for q = 3/2,7/5,4/3 (Fig. 8c) and for a time lag At = 40 days for 
q = 7/5, 4/3, 5/4 (Fig. 8d). As expected the tails of the distribution functions 
approach a Gaussian type when q is approaching 1. For completeness, the 
corresponding cases of the distribution of price increments are shown and 
briefiy discussed in the Appendix. 

In so doing the probability distributions of the normalized log returns for 
the different values of the time lag At = 1, 5, 10, 15, 20, 25, 30, 35, 40 days can 
be shown in Fig. 2 together with the lines representing the best fit to the 
Tsallis type of distribution function. In Table 3 the statistical parameters 
related to the Tsallis type of distribution function are summarized, including 
a criterion for the goodness of the fit, i.e. the Kolmogorov-Smirnov distance 
dKs, which is defined as the maximum distance between the cumulative prob- 
abihty distributions of the data and the fitting lines. Note that the kurtosis 
(see Table 3) for the TsaUis type of distribution function 

where Kl = 3 for a Gaussian process, is positive for all values of q < 7/5 
as expected, since its positiveness is directly related to the occurrence of 
intermittency [6]. Moreover, the hmit g < 7/5 also imphes that the second 

^Sattin formula [23] might also be tested in future work. 
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moment of the Tsallis type distribution function will always remain finite, 
as necessarily due in the type of phenomena hereby studied. Furthermore, 
if we assume that the Kolmogorov log-normal model of turbulence [47] is 
applicable and let At^^ be the scale at which the whole partial distribution 
function becomes Gaussian, then the kurtosis Kr should scale as 

K,^K.{^f. (12) 

Therefore 

,~7(At/At,r (13) 

In order to obtain an estimate for At/,, we observe that the turbulence 
model, Eq. (4), fits well the normalized log returns for At = 120 days and q = 
1.01 (Fig. 9a). The a-parameter {a — 0.74) in this case plays an important 
role in controlling the tails such that the Tsallis type distribution function 
for negative values of Z fits the data whose probability distribution function 
still deviates from Gaussian. In fact, further increasing the time lag to the 
value At = 200 days leads to a complete coincidence between the distribution 
functions in the Tsallis and Gaussian forms for the presently investigated data 
(Fig. 9b). Corresponding parameter values are also listed in Table 3. This 
short observation convincingly indicates where the transition occurs between 
the small time scale model of nonextensive, intermittent process and the large 
scale Gaussian extensive homogeneous fluctuation picture [6, 15] and refine 
the estimate of the Gaussian range in Figs. 1-2. 

In Fig. 10 the Tsallis parameter q is shown as a function of the rescaled 
time lags At/At j^, where Ati is the integral scale, the scale at which the 
whole probability distribution function converges to Gaussian. The crosses) 
represent the q values for which the best fit to the S&P500 data (Fig. 2) 
is obtained with Eq. (4). With this the value of the integral scale At^, we 
find the value of the exponent S — 0.39 as the one for which the Eq. (13) 
fits best the g-values. The exponent value S = 0.39 also allows to fit well the 
power law dependence (Eqs. (11) and (12)) of the rescaled kurtosis K^/ Kl 
as shown in the inset of Fig. 10. 

Note that in the framework of the Kolmogorov log- normal model [47, 20], 
5 — Ak/9, where k is called the intermittency exponent. Therefore, we find 
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K — 0.88 for the intermittency exponent of normalized log returns of the 
S&P500 daily closing price in the time interval of interest. This value of k is 
higher than the value of the intermittency exponent k — 0.25 for turbulence 
recently obtained from experimental atmospheric data [48]. Early estimates 
have varied from 0.18 to 0.85 using different experimental techniques [49, 50, 
51]. Large values of the intermittency exponent, ranging from 0.2 to 0.8, have 
been reported in studies of multiparticle production [52]. It was found that 
the range of intermittency exponent values depend on the number of cascades; 
the smaller the number of stages of the multiplicative cascade the smaller n, 
and conversely [Fig. 2b in [52]]. In analogy with such findings, a value of 
K = 0.88 can be considered to be related to a high number of cascades in a 
multiplicative process, leading to the observed partial distribution functions 
of the normahzed log returns of the S&P500 index. 

One can explore the Tsallis type of the probability distribution function 
Eq.(4) in two limits. For small values of normalized log returns Z the prob- 
ability distribution function converges to the form 

Therefore the Tsallis type distribution function converges to a Gaussian, i.e. 
a — > 1, for small values of the normalized log returns, for any At investi- 
gated hereby (see Figs. 1-2). It is also of interest to check the probability of 
return to the origin pAt (-2^ = 0) (Table 3). There is a slight difference between 
the values of the probabihty of "return to the origin" for the data and the 
one obtained from Eq.(4) pAt{Z = 0) = i/Zg. This difference is decreasing 
with increasing At and completely disappears in the Gaussian limit q — > 1, 
a — > 1. 

In the limit of large values of normalized log returns Z, the Tsallis type 
distribution converges to a power law 




q-l- 



2a -{q- 1) 

Studying the Tsallis type of distribution function one can obtain from 
Eq. (4) an expression for the width of the Tsallis type of probability distribu- 
tion function, 2a^ = {2a — {q — l))/{2aCPo{q — 1)). In the hmit of a — >• 1 
the width of the Tsallis type distribution 2(7^ = (3 — q)/2C/3o{q — 1), i.e. 
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~ 2/(C/9o). It is obvious that for large time lags 2f7^, tends to diverge [24], 
like ~ (At)^/(^~'^); this can be easily verified on a log- log plot (not shown). 

In limit of q — > 1 the Tsallis type distribution function converges to 
Gaussian (as seen in Fig. 9b). The values of the parameters q, a, C^q, that 
best fit the data using Eq.(4), and 2(7^ are plotted as a function of the time 
lag in Fig. 11. 



5 Fokker-Planck approach 

On the other hand, the evolution of a time dependent probability distribu- 
tion function is usually described within the Fokker-Planck approach. This 

method provides some further information on the correlations present in the 
time series and it begins with the joint PDF's, that depend on A/" variables, 
i.e. p^{Zi, Ati] At_A/)- We started to address this issue by determin- 
ing the joint PDF for M — 2, i.e. p{Z2, At2\ Axi, Ati). The symmetrically 
tilted character of the joint PDF contour levels (Fig. 12) around an inertia 
axis with slope 1/2 points out to some statistical dependence, i.e. a correla- 
tion, between the normalized log returns Z(t, At) of the daily closing price 
signal of S&P500. A lack of correlations would put the inertia axis on the 
main diagonal (Fig. 12). 

The conditional probability function is 

p(Z.^„ AU^^\Z,, AU) = -^^-—-^ (16) 

for i — 1,...,A/'— 1. For any At2 < Ati < Ati, the Chapman-Kolmogorov 
equation is a necessary condition of a Markov process, one without memory 
but governed by probabilistic conditions 



p(Z2, At2|^i, Ah) = I d{Zi)p{Z2, At2\Zi, Ati)p{Axi, AU\Zi, Ah). (17) 

The Chapman-Kolmogorov equation when formulated in dif ferential 
form yields a master equation, which can take the form of a Fokker-Planck 
equation [43]. Let r = log2{200 / At) , 



d 

dr 



p{Z,t) 



p{Z,t) 



(18) 
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in terms of a drift D'^^\Z,t) and a diffusion coefficient D^'^\Z,t) (tlius values 
of r represent Ati, i = I, ■■■)■ 

The coefficient functional dependence can be estimated directly from 
the moments M^''^ (known as Kramers-Moyal coefficients) of the conditional 
probability distributions: 

= -^J dZ\z' - Zfp{Z\T + At|Z, t) (19) 

L'('=)(Z,r) = ^limM^'^) (20) 

for Ar — > 0. According to Fig. 13a the drift coefficient D^^^ and the 
diffusion coefficients I?^^) is well represented (Fig. 13b) by a parabola 

D(2)(^) = 0.26^2 - 0.005Z + 0.02 (21) 

in the interval Z G [—0.175, 0.225], - noticing that it is smaller than the one 
presented in Fig. 2. 

It may be worthwhile to recall that the observed quadratic dependence 
of the diffusion term D^^) is essential for the logarithmic scahng of the inter- 
mittency parameter in studies on turbulence. 

Finally, the Fokker-Planck equation for the distribution function is known 
to be equivalent to a Langevin equation for the variable, i.e. Z here, (within 
the Ito interpretation [7, 43, 44, 45, 46]) 

^Z(r) = L>«(Z(r),r)+r^(r)Vi?(2)(Z(r),r), (22) 

where rjir) is a fluctuating 5-correlated force with Gaussian statistics, i.e. < 
77(r) r]{T')> = 25{t-t'). 

Thus the Fokker-Planck approach provides the evolution process of PDF's 
from small time lags to larger ones. The fact that the drift coefficient is 
approximately equal to zero, therefore indicating that there is no correlation 
between the probability distribution functions for different time lags, is well 
related to the Gaussian character of the distribution function for such small 
values of the normalized log returns Z G [—3, 3]. D^^^ ^ further implies that 
there is almost no "restoring force", i.e. 7 ^ in Eq. (6), while the quadratic 
dependence of D'^'^^ in Z is obviously like an autocorrelation function for a 
diffusion process. 
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6 Conclusion 



In summary, we have presented a method that provides the evolution pro- 
cess of probabiUty distribution functions (over twenty years) of one financial 
index, i.e. the S&P500. We have studied the evolution process of the tails 
that are outside the central (Gaussian) regime at small returns, thereby fa- 
cilitating the understanding of the evolution of these distribution functions 
in a Fokker-Planck framework. Beck turbulence model can be well applied 
to describe the volatility (of normalized log-returns) distributions assuming 
a x^-distribution for the " local" volatility. An open question in nonextensive 
thermostatistics studies is often raised about the meaning, value and behav- 
ior of the non extensive exponent, or Tsallis parameter q. The inter mittency 
exponent is found to be related to the scaling exponent of the PDF moments 
in the framework of Kolmogorov log-normal model, thereby giving weight to 
the model and the statistical approach. The intermittency exponent large 
value points to a large number of cascades in the turbulent process. Its range 
has been found to extend up to ca. 200 days. One may still wonder on the q- 
value itself. In other works, this value is related, e.g. to the upper and lower 
bounds of the multifractal dimension [17], in other words to the bounds of 
the a values in multifractal studies [40]. It may also be related to the value 
of the fractional derivative, say in a non-linear Fokker-Planck equation ap- 
proach [53]. This should be some interesting work to pursue, - again with 
some warning concerning the possible error bars on the generalized fractal 
dimension in multifractal studies [54]. 

We have also presented the turbulence-like dynamics through the Fokker- 
Planck and the Langevin equations. We have (unexpectedly) found that, in 
the treated case, there is almost no "restoring force", i.e. (7 in the 
Langevin equation). A comparison is made between normalized log-returns 
and mere price increments. We have examined the corresponding cases of 
the distribution of price increments with other possible definitions. It was 
found that the definition (through a normalized log return rather than a 
mere price difference) is very relevant for obtaining nice fits. This has been 
also observed in a work by Karth and Peinke [55] on related matter. This 
warning might also shine some light on the possible origin of the controversy 
[4, 5] concerning the relationship (or not) between the fat tails caused by 
some dynamical hierarchical cascade process of volatility correlations. 

These points not withstanding, we have related a financial market behav- 
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ior to Tsallis non extensive thermodynamics approach, i.e. more precisely to 
a turbulence-hke process, - as financial market and indices were often claimed 
to be seen [2, 12, 13]. Finally, it seems that we have thoroughly answered the 
often raised question "why to look at the tails of a probability distribution 
function? and what does that lead to?" . 

7 Appendix 

We have also searched for describing the partial distribution function of the 
(raw) increments of daily closing price signal of the S&P500 with the TsaUis 
type distribution function. We have applied Eq.(4) for x = Ay(t,At) — 
y{t + At) — y{t). We have tested the Tsallis type distribution function for 
the increments of At = 1 day for fixed q = 1.45 {Cf3o = 0.23) and varying 
a = 0.5, 0.6, 0.7 (Fig. 14a). Applying the same set of parameters and to price 
increments for Ay = 40 days, leads to a pretty bad fit (Figs. 14b). Decreasing 
the value of q would not have produce better results since the Tsallis type 
distribution function would have been bounded within smaller range around 
Ay = values. A test for fixed a = 0.5 {C(3o = 0.23) and varying q = 
1.45, 1, 30, 1.15 for At = 1 day is next shown in Fig. 14c. Again the same set 
of parameters is applied to price increments for Ay = 40 days and lead to a 
pretty bad fit (Figs. 14d). These results may be somewhat expected because 
the Tsallis type distribution function represents a mathematical construction 
that is designed for normalized variables, i.e. a variable changing within a 
limited range. To take into account a double peak like structure (e.g. for 
large time lags, see Fig. 14) remains an open question. 
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Figure Captions 



Figure 1 - Probability distribution function of normalized log returns Z{t, At) 
of daily closing price value signal of S&P500 between Jan. 01, 1980 and Dec. 
31, 1999, for At = 1 day (symbols). Normalized log returns are calculated as 
Z{t,At) = iyit)- < y >At)/f^At, where y{t) = ln{y{t + At) / y{t)) and ciAt is 
the standard deviation of y{t) for time lag At. The dashed line represents a 
Gaussian distribution. Inset: Power law fit (solid hue) of the negative (-3.1) 
and positive (-3.1) slope of the distribution outside the Gaussian regime 

Figure 2 - Probability distribution function (PDF) p^t{Z) of normalized 
log returns of the daily closing price value signal of S&P500 (symbols) and 
the Tsallis type distribution function (lines) for different values of At — 
1, 5, 10, 15, 20, 25, 30, 35, 40 days. The PDF (symbols and curves) for each At 
are displaced by 10 with respect to the previous one; the curve for At = 1 day 
is unmoved. The large circles mark the ends of the interval in which the 
distribution is like a gaussian distribution. The values of the slopes of the 
positive and negative tails of the distributions outside the gaussian range 
Z e [—3, 3] are listed in Table 1. The values of the parameters for the Tsallis 
type distribution function for each At are summarized in Table 3 

Figure 3 - DFA function F{n) plotted as a function of the the box size 
n of the integrated daily closing price value signal of S&P500 between Jan. 
01, 1980 and Dec. 31, 1999. Brownian-like fluctuations with 1 -|- H^fa — 
1.52 ± 0.01 are obtained for all possible time scales 

Figure 4 - DFA function F{n) plotted as a function of the the box size n 
of the integrated normahzed log returns Z(t, At) of daily closing price value 
signal of S&P500 between Jan. 01, 1980 and Dec. 31, 1999 for different time 
lags At = 1, 5, 10, 15, 20, 25, 30, 35, 40 days. Values of the scaling exponents 
Hdfa for the various DFA functions are summarized in Table 2. Inset: 
Power law functional dependence of the value of the cross over box-size Ux 
as a function of time lag At 

Figure 5 - Power spectrum S{f) of the daily closing price value signal of 
S&P500 between Jan. 01, 1980 and Dec. 31, 1999. A scale break at around 
/ = 1/250 day^^ separates two scaling regions. Inset: Scaling of the power 
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spectrum of the daily closing price signal of S&P500 as a white noise signal 
with n ^ 

Figure 6 - Power spectrum S{f) of the normalized log returns Z(t, At) 
of daily closing price value signal of S&P500 between Jan. 01, 1980 and 
Dec. 31, 1999 for different time lags At = 1, 5, 20,35,40 days. Each curve is 
displaced by 10^^ with respect to the previous one; the power spectrum of 
the normalized log returns for At — 1 day is not displaced. The dashed hne 
from / = 1/70 days~^ to / = 1/2 days~^ has a slope /i = 1.86, corresponding 
to the Hdpa exponent. The horizontal dashed line from / = 1/10^'^ days~^ 
to / = 1/128 days~^ corresponds to what should be expected for white noise 
and is in agreement with the scaling of the DFA function for the same data 
in Fig. 4 

Figure 7 - Probability density f'AtiP) of the local volatility j3 (Eq.(lO)) in 
terms of standard deviation of the normalized log returns Z(t. At) of S&P500 
in non-overlapping windows with size m=32 days for different time lags (sym- 
bols) (a-i) At =1, 5, 10, 15, 20, 25, 30, 35, 40 days. Lines: x^-distribution 
as given by Eq. (7). 

Figure 8 - Probability distribution functions of the normalized log returns 
of daily closing price signal of S&P500 (symbols) for (a) At = 1 day and 
fixed q — 1.39. The TsaUis type distribution functions (Eq.(4)) obtained for 
various values of the parameter a — 1.0,0.9,0.8, dashed, solid, dash-dotted 
line, respectively; (b) same as (a) but for At = 40 days and g = 1.22; (c) for 
At = 1 day and fixed a = 1.0 for various values of g = 3/2, 7/5, 4/3, dashed, 
solid, dash-dotted line, respectively; (d) same as (c) but for At = 40 days 
and g = 7/5,4/3,5/4 

Figure 9 - Partial distribution function of normalized log return of daily 
closing price of S&P500 for a large time lag, i.e. (a) At = 120 days and 
(b) At = 200 days. The solid line marks the best fit with a TsaUis type 
distribution function, Eq. (4), while the Gaussian distribution function is 
drawn with a dashed line 

Figure 10 - The functional dependence of the Tsallis q parameter on the 
rescaled time lag At/Ati for AtL = 200 days and 6 = 0.39 (see Eq. (13)) 
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(line); the symbols represent the values of the q parameter listed in Table 3 
and used to plot the fitting lines in Figs. 2 and 9. Inset : Scaling properties 
of the rescaled kurtosis K^/Kl, where = 3 is the kurtosis for a Gaussian 
process, as a function of the rescaled time lag Ai/Ai^ satisfying Eq. (11) 
(open symbols) and Eq. (12) (full symbols) 

Figure 11 - Characteristic parameters of Tsallis type distribution function 
as defined in [27] : TsaUis ^-parameter (crosses), a (squares), constant C/3o 
used in the fit (open circles), the width of the Tsallis type distribution 2o"^ = 
{2a — {q — l))/(2«C/3o(q' — 1)) from Eq.(4) (triangles) (rescaled by a factor 
of 1/6), asymptotic behavior of 2(T^ a; 2/(C/3o) for a — > 1 (full circles) 
(rescaled by a factor of 1/6) 

Figure 12 - Typical contour plots of the joint probability density function 
p{Z2, At2; Zi, Ati) of daily closing price of S&P500 for the period of interest 
Jan. 01, 1980 and Dec. 31, 1999. Dashed lines have a slope +1 and emphasize 
the correlations between probability density functions for (a) At2 = 1 day 
and Ati — 5 days and (b) At2 — 5 days and A^i = 10 days. Contour 
levels correspond to logwp{Z2, A.t2, Zi, Ati) = —1.0, —1.5, —2.0, —2.5, —3.0 
from center to border 

Figure 13 - Kramers-Moyal drift and diffusion coefficients (a) D'^^^ and 

(b) D^^^ as a function of normalized log returns Z for daily closing price of 
S&P500 ; ^ 0.26^2 - 0.0005Z + 0.02 

Figure 14 - Probabifity distribution functions of the daily closing price 
increments Ay{t) = y{t + At) — y{t) of S&P500 (symbols). The Tsallis type 
distribution functions (Eq.(4)) obtained (a) for At — 1 day and fixed q — 1A5 
{CPo = 0.23) and for various values of the parameter a = 0.5, 0.6, 0.7, dashed, 
solid, dash-dotted line, respectively; (b) the same as (a) but for At = 40 days; 

(c) for At = 1 day and fixed a = 0.5 (C/3o = 0.23) and for various values 
of g = 1.45, 1.30, 1.15, dashed, solid, dash-dotted line, respectively; (d) the 
same as (b) but for At — 40 days 
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Table 1: Slopes of the positive and negative tails (second and third column) 
of the distributions outside the Gaussian range Z{t,At) e [—3,3]; values of 
characteristic Tsallis function parameters are given 



At 


positive 


negative 


2q 


1) 


a 


1 


3.1 


3.1 


2.78 


2.564 


0.92 


5 


3.1 


3.0 


2.72 


2.778 


0.90 


10 


3.0 


2.7 


2.68 


2.941 


0.88 


15 


3.0 


2.6 


2.66 


3.030 


0.86 


20 


2.9 


2.5 


2.64 


3.125 


0.85 


25 


2.9 


2.5 


2.62 


3.226 


0.83 


30 


2.8 


2.3 


2.58 


3.448 


0.80 


35 


2.7 


2.3 


2.50 


4.000 


0.78 


40 


2.5 


2.2 


2.44 


4.545 


0.76 



Table 2: Values of the scaling exponent from the DFA analysis of 
normalized log returns Z for different values of the time lag At = 
5, 10, 15, 20, 25, 30, 35, 40 days, and crossover "box size" Ux 



At 


1 + HoFA 


nx 


5 


1.27±0.04 


19 


10 


1.37±0.02 


35 


15 


1.39±0.02 


49 


20 


1.38±0.02 


70 


25 


1.38±0.02 


91 


30 


1.41±0.01 


108 


35 


1.43±0.01 


117 


40 


1.43±0.01 


128 
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Table 3: Values of the parameters characterizing the S&P500 daily closing 
price data between Jan. 1, 1980 and Dec. 31, 1999 in the nonextensive 
thermostatisitics approach. For the definition of the Kolmogorov-Smirnov 
distance dxs see the text 



At 
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a 


CPo 


PAtiZ = 0) 


PAtiZ = 0) 


Kr 


dxs 














data 


Eq.(4) 
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.39 


0. 


,92 


0.65 


0.505 


0.611 


49.800 


0.072 


5 


1 


.36 


0. 


,90 


0.62 


0.447 


0.600 


13.800 


0.100 


10 


1 


.34 


0. 


,88 


0.60 


0.462 


0.592 


9.800 


0.091 


15 


1 


.33 


0. 


,86 


0.58 


0.472 


0.582 


8.657 


0.085 


20 


1 


.32 


0. 


,85 


0.56 


0.459 


0.572 


7.800 


0.085 


25 


1 


.31 


0. 


,83 


0.54 


0.447 


0.560 


7.133 


0.087 


30 


1 


.29 


0. 


,80 


0.52 


0.443 


0.549 


6.164 


0.088 


35 


1 


.25 


0. 


,78 


0.50 


0.432 


0.538 


5.000 


0.087 


40 


1 


.22 


0. 


,76 


0.48 


0.445 


0.525 


4.467 


0.077 


120 


1 


.01 


0, 


,74 


0.39 


0.431 


0.467 


3.031 


0.052 


200 


1, 


.01 


1. 


,00 


0.26 


0.398 


0.406 


3.031 


0.040 
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Figure 1: Probability distribution function of normalized log returns Z{t, At) 
of daily closing price value signal of S&P500 between Jan. 01, 1980 and Dec. 
31, 1999, for At = 1 day (symbols). Normalized log returns are calculated as 
Z{t,At) = {y{t)- < y >At)/crAt, where y{t) = ln{y{t + At)/y{t)) and is 
the standard deviation of y{t) for time lag At. The dashed line represents a 
Gaussian distribution. Inset: Power law fit (solid line) of the negative (-3.1) 
and positive (-3.1) slope of the distribution outside the Gaussian regime 
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Figure 2: Probability distribution function (PDF) pAt{Z) of normalized 
log returns of the daily closing price value signal of S&P500 (symbols) 
and the Tsallis type distribution function (lines) for different values of 
At = 1, 5, 10, 15, 20, 25, 30, 35, 40 days. The PDF (symbols and curves) for 
each At are displaced by 10 with respect to the previous one; the curve for 
At = 1 day is unmoved. The large circles mark the ends of the interval 
in which the distribution is like a gaussian distribution. The values of the 
slopes of the positive and negative tails of the distributions outside the gaus- 
sian range Z G [—3, 3] are listed in Table 1. The values of the parameters for 
the Tsallis type distribution function for each At are summarized in Table 3 
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Figure 3: DFA function F{n) plotted as a function of the the box size n of the 
integrated daily closing price value signal of S&P500 between Jan. 01, 1980 
and Dec. 31, 1999. Brownian-like fluctuations with 1 + H^tpj^ = 1.52 ± 0.01 
are obtained for all possible time scales 
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Figure 4: DFA function F{n) plotted as a function of the the box size n of the 
integrated normahzed log returns Z{t, At) of daily closing price value signal 
of S&P500 between Jan. 01, 1980 and Dec. 31, 1999 for different time lags 
At = 1, 5, 10, 15, 20, 25, 30, 35, 40 days. Values of the scaling exponents Hdfa 
for the various DFA functions are summarized in Table 2. Inset: Power law 
functional dependence of the value of the cross over box-size as a function 
of time lag At 



30 




1/10000 1/1000 1/100 1/10 1 



f (1/day) 

Figure 5: Power spectrum S{f) of the daily closing price value signal of 
S&P500 between Jan. 01, 1980 and Dec. 31, 1999. A scale break at around 
/ = 1/250 day~^ separates two scaling regions. Inset: Scaling of the power 
spectrum of the daily closing price signal of S&PSOO as a white noise signal 
with /i ^ 
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Figure 6: Power spectrum S{f) of the normalized log returns Z{t, At) of 
daily closing price value signal of S&P500 between Jan. 01, 1980 and Dec. 
31, 1999 for different time lags At = 1,5,20,35,40 days. Each curve is 
displaced by 10^^ with respect to the previous one; the power spectrum of 
the normalized log returns for At = 1 day is not displaced. The dashed line 
from / = 1/70 days~^ to / = 1/2 days~^ has a slope fi = 1.86, corresponding 
to the Hijpji exponent. The horizontal dashed line from / = 1/10"^ days~^ 
to / = 1/128 days~^ corresponds to what should be expected for white noise 
and is in agreement with the scaling of the DFA function for the same data 
in Fig. 4 
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Figure 7: Probability density fAtiP) of the local volatility P (Eq.(lO)) in 
terms of standard deviation of the normalized log returns Z{t, At) of S&P500 
in non-overlapping windows with size m=32 days for different time lags (sym- 
bols) (a-i) At =1, 5, 10, 15, 20, 25, 30, 35, 40 days. Lines: x^-distribution 
as given by Eq. (7) 
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Figure 8: Probability distribution functions of the normalized log returns 
of daily closing price signal of S&P500 (symbols) for (a) At = 1 day and 
fixed q = 1.39. The Tsallis type distribution functions (Eq.(4)) obtained for 
various values of the parameter a = 1.0, 0.9, 0.8, dashed, solid, dash-dotted 
line, respectively; (b) same as (a) but for At = 40 days and g = 1.22; (c) for 
At = 1 day and fixed a = 1.0 for various values of g = 3/2, 7/5, 4/3, dashed, 
solid, dash-dotted line, respectively; (d) same as (c) but for At = 40 days 
and q = 7/5,4/3,5/4 
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Figure 9: Partial distribution function of normalized log return of daily clos- 
ing price of S&P500 for a large time lag, i.e. (a) At = 120 days and (b) 
At = 200 days. The solid line marks the best fit with a Tsallis type distri- 
bution function, Eq. (4), while the Gaussian distribution function is drawn 
with a dashed line 
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Figure 10: The functional dependence of the Tsalhs q parameter on the 
rescaled time lag At/AtL for Ati = 200 days and 6 = 0.39 (see Eq. (13)) 
(line); the symbols represent the values of the q parameter listed in Table 3 
and used to plot the fitting lines in Figs. 2 and 9. Inset : Scaling properties 
of the rescaled kurtosis K^/Kl, where = 3 is the kurtosis for a Gaussian 
process, as a function of the rescaled time lag At/Ati satisfying Eq. (11) 
(open symbols) and Eq. (12) (full symbols) 
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Figure 11: Characteristic parameters of Tsallis type distribution function 
as defined in [27] : Tsallis g-parameter (crosses), a (squares), constant C^q 
used in the fit (open circles), the width of the Tsallis type distribution 2cr^ = 
(2a — (g — l))/(2Q;C/3o(g — 1)) from Eq.(4) (triangles) (rescaled by a factor 
of 1/6), asymptotic behavior of 2(T^ ^ 2/{C(3q) for a — > 1 (full circles) 
(rescaled by a factor of 1/6) 
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Figure 12: Typical contour plots of the joint probability density function 
^(^2, At2; Ati) of daily closing price of S&P500 for the period of interest 
Jan. 01, 1980 and Dec. 31, 1999. Dashed lines have a slope +1 and emphasize 
the correlations between probability density functions for (a) At2 = 1 day 
and Ati = 5 days and (b) At2 = 5 days and Ati = 10 days. Contour 
levels correspond to logiop{Z2, At2; Zi, Ati) = —1.0, —1.5, —2.0, —2.5, —3.0 
from center to border 
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Figure 13: Kramers- Moyal drift and diffusion coefficients (a) D^^^ and (b) 
function of normalized log returns Z for daily closing price of 
S&P500 ; ^(2) = 0.26^2 - 0.0005Z + 0.02 
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Figure 14: Probability distribution functions of the daily closing price in- 
crements A?/(t) = y{t + At) — y{t) of S&P500 (symbols). The Tsallis type 
distribution functions (Eq.(4)) obtained (a) for At = 1 day and fixed q = 1.45 
{CPo = 0.23) and for various values of the parameter a = 0.5, 0.6, 0.7, dashed, 
solid, dash-dotted line, respectively; (b) the same as (a) but for At = 40 days; 
(c) for At = 1 day and fixed a = 0.5 {CPo = 0.23) and for various values 
of g = 1.45,1.30,1.15, dashed, solid, dash-dotted line, respectively; (d) the 
same as (b) but for At = 40 days 



40 



